Schrodinger cat states prepared by Bloch oscillation in a spin-dependent optical lattice 

B. J. Wu and J. M. Zhang 

Beijing National Laboratory for Condensed Matter Physics, 

Institute of Physics, Chinese Academy of Sciences, Beijing 100080, China 

We propose to use Bloch osciUation of ultra-cold atoms in a spin-dependent optical lattice to 
prepare schrodinger cat states. Depending on its internal state, an atom feels different periodic 
potentials and thus has different energy band structures for its center-of-mass motion. Consequently, 
under the same gravity force, the wave packets associated with different internal states perform Bloch 
oscillation of different amplitudes in space and in particular they can be macroscopically displaced 
with respect to each other. In this way, a cat state can be prepared. 

PACS numbers: 03.75.Lm, 03.65.Ud, 03.65.Sq 
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Bloch oscillation is a peculiar response of a particle 
in a periodic potential to a weak external force [l|, |2]- 
Under the drag of a constant force F, the wave packet 
of the particle oscillates back and forth periodically in 
space without accelerating indefinitely in the direction 
of the force. The weirdness is that a DC bias generates 
an alternating current. The underlying reason is that as 
long as the force is weak enough, interband transitions 
are prohibited by the gaps between energy bands. Con- 
fined in a specific energy band by such a mechanism, the 
motion of the particle is captured to a good extent by 
the semi-classical equations [3!] 
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where r and q are the center-of-mass and wave vector 
of the particle, respectively, and En{q) is the dispersion 
relation in the n-th band. In the g-space, the picture is 
that the particle traverses the Brillouin zone, which is 
of the topology of a circle, repeatedly at a constant rate. 
From these equations, one solves readily the displacement 
of the wave packet as 



1 



Ar(i) = - [En{qo + Ft/h) - £;„(qo)] 
r 



(2) 



where q^ is the initial value of q. We see that Ar is a peri- 
odic function of time with the period T — 27Th/aF — this 
states the Bloch oscillation in the semi-classical theory. 
Here a is the period of the periodic potential and 27r/a 
is the size of the Brillouin zone. Note that T is indepen- 
dent of the detailed structure of the periodic potential 
but depends only on its period. 

Suppose initially the particle is at the bottom of the 
lowest energy band (n = 0) with go = 0- The maximum 
displacement Ar^ is reached at t = T/2, when the wave 
vector q — n/a arrives at the edge of the first Brillouin 
zone. Afterwards, the velocity of the particle dr/dt re- 
verses. The value of Ar™ is simply 
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Ar,„ = — , 



(3) 



where B = Eo{TT/a) — Eo{0) is the band width of the 
lowest band. Simple as it looks, this equation has an 
important implication in our work below. 



So far, Bloch oscillation has been observed in a variety 
of systems, such as semiconductor super lattices [J, [5|, 
cold atoms in optical lattices |6l-[llj. and photonic lat- 
tices |12| - [l4| . On the application side, it has found use 
in microwave generation [J, |5|, llSl ] , precision force mea- 
surements jSMllj , and coherent transport of matter waves 

[IMl. 

In this paper, we propose that Bloch oscillation can 
also be used to prepare Schrodinger cat states [l9l - [23| | 
in a spin-dependent optical lattice. The idea is actually 
very simple. In a spin-dependent optical lattice, atoms 
in different internal states see different potentials (e.g. 
of different strengths) . This non-trivial fact means that 
they also have different energy band structures for their 
center-of-mass motion, which in turn means they will 
have different Bloch oscillation modes under the same 
force (see Eq. ([U). In particular, their maximum dis- 
placements will be different according to ([3]) since the B's 
may differ, and this implies that the wave packets cor- 
responding to different internal states will be displaced 
with respect to each other. That is, a cat state can be 
prepared. 

First of all we need a spin-dependent optical lattice. 
As proposed in [2J, [23 and realized in [2^ , such an opti- 
cal lattice can be constructed by interfering two counter- 
propagating laser beams linearly polarized but with an 
angle 9 between the polarization vectors. The resulting 
standing light field can be decomposed into a cr+ and a 
(T_ polarized one with intensities /+ = Im cos^{kx + 6/2) 
and /_ = Im cos^{kx — 6/2), respectively. Here k = 2tt/X 
is the wave vector of the laser beams. Such a decomposi- 
tion is helpful since the dipole potential for an atom in a 
state \F,mp) (the quantization axis of the atom is along 
the optical lattice) is simply the sum of the contributions 
of the two components. Below we will use the same sys- 
tem as in [2g|. That is, we choose ^^Rb as the atom, and 
1) = \F = 2,mF = -2) and |0) = \F = l,mF = -1) as 
the two atomic internal states. The spin dependence of 
the dipole potential is realized by choosing the laser fre- 
quency to resolve the fine structure of the Rubidium D 
line. Specifically, as in [26|, by tuning the wave length of 
the optical lattice laser to A = 785 nm, the dipole poten- 
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FIG. 1: (Color online) (a) The dipole potentials Vi{x) (solid line) and Vo{x) (dashed line) for an atom in the internal states 
|1) and |0), respectively. The parameters are Vm/Er = 5 and 9 = ix 12. Here Er = K^k^ /2m is the recoil energy of the atom, 
(b) The lowest three energy bands for an atom in the |1) (solid lines) and |0) (dashed lines) internal states. 



tials for an atom in the |1) and |0) states are respectively, 



Vi{x\9) = Vmcos [kx- 



2/' 



Vo{x; 9) = -V„, cos^ (kx+-\+ -V„, cos^ ( '^'^ ~ 2 ) ' 

where Vm oc Im- The point is that if 7^ or tt, Vi{x) 
and Vo{x) are shifted relative to each other, and more 
importantly, have different amplitudes. The very latter 
effect results in different band structures as we see in 
Fig. [TJ There the potentials Vi^ and the corresponding 
energy bands for the two internal states are depicted. 
The parameters chosen are Vm/Er = 5 and = tt/2, 
where E,. = h^k"^ /2m = 27rh x 3.72 kHz is the recoil 
energy of the atom. Note that the width of the lowest 
band for the |1) state is 2Trh x 0.983 kHz, while that for 
the |0) state is 2TTh x 1.925 kHz. The two differ almost 
by a factor of 2. 

Now our scheme to generate a Schrodinger cat goes like 
this. Suppose initially the angle 6 = (for this value of 
9, Vi = Vq) and the atom is in the |1) state. As for its 
external state, it is assumed to be 
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dqf{q)(j)o{q)- 
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Here (j)o {q) is the Bloch state in the lowest band with wave 
vector q. The weight function f{q) is localized around 
q = but otherwise unspecified. This condition is easily 
satisfied as long as the spatial size of the wave packet ^^ 
is much larger than the lattice constant a = A/2. Ac- 
tually, the condensate wave function should satisfy this 
condition if the condensate is loaded adiabatically from 
a magnetic trap into the lattice as is usually done in cold 
atom experiments. It is checked that the results pre- 
sented in the following are barely affected with different 
choices of /(g) , as long as the localization condition is sat- 
isfied. This fact is consistent with the semi-classical the- 
ory in which the details of the wave packets are irrelevant. 



Specifically, in the simulations to be presented, f{q) is of 
the form f{q) oc exp{~q^/w^) with wa/n = 0.1 ^ 1. 
This value of w corresponds to a wave packet with a size 
on the order of 10a. 

Then at some moment, by using a microwave pulse 
we can prepare the internal state of the atom into an 
arbitrary superposition of the |0) and |1) states. Let it 
be a|0) -t- /3|1) with \a\'^ + |/3|^ = 1. Subsequently, 9 is 
adjusted to 7r/2 (suddenly or smoothly, it does not mat- 
ter; but in our simulation we take the sudden scenario), 
and the lattice is tilted by ip with respect to the hor- 
izontal plane. The Bloch oscillation then starts. At a 
later time i, the wave function of the atom is of the form 
^(i) =a|0)*o(a;,i) +/3|l)*i(a:,i)- The evolution of the 
external wave function ^j(a;,i) is given by (j = 0, 1) 






2ni dx^ 



+ Vj{x)^Fx] * 



(5) 



with the force F = nig sin ip and the initial condition 
*,(i = 0) = *,. 

Although the semi-classical theory above gives us an 
overall idea of the motion of the wave packets ^j(x,t), 
here we shall solve Eq. ([5|) numerically. Snapshots of 
^j(a:,i) are shown in Fig. [5] As expected, in the inter- 
val < i < T/2, the two wave packets move rightward 
and gradually depart. At the turning point t = T/2 (see 
Fig.HJl), the distance between the two reaches the max- 
imum. Remarkably, at this point, the two wave pack- 
ets are well separated — the separation between them is 
about 50a, which is much larger than their sizes (~ 10a). 
Note that the wave packets themselves are large enough 
to deserve the name macroscopic. Therefore we obtain a 
desired Schrodinger cat state 

*(i == T/2) = a|l)*i(i = T/2) + /3|0)*o(i = T/2), (6) 

for which the internal and external states of the atom 
are entangled. The point is that the latter is in two 
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(f) t/T= 1 
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FIG. 2: (Color online) Time evolution of the wave packets 'i/i{x,t) (blue lines) and 'i/o{x,t) (green lines) in a Bloch oscillation 
cycle. The parameters are the same as in Fig. [T]and Fa/Er — 0.005. Initially, 'i/ifi — ^i as defined in Q. Note that the wave 
packets are not normalized to unity. 



macroscopically displaced macroscopic states, which cor- 
respond to a "live" and a "dead" cat, respectively. The 
coherence between the two packets can be checked by 
applying a microwave pulse to achieve a rotation in the 
{|1), |0)} space, and then turning off the lattice and ob- 
serving the momentum distribution of the atom in the 
|1) state by the absorption imaging method I26l]. Note 
that ^1,0 (r/2) are both peaked around q = ±7r/a in the 
momentum space and thus the interference pattern will 
primarily consist of two peaks whose amplitudes depend 
on the rotation as well as a and /3. 

Afterwards, the two wa ve p ackets move backwards 
and return approximately [27| to their original states 
at t = T (see Fig. [5R That is, *i,o(r) ~ *» up to 
some global phases [2g]. The global wave function is 
then *(T) ~ (a|0) -h /^e^'^ll))*,, where x is the differ- 
ence of the phases ^i,o accumulated in a cycle. Thus 
the atom completes a Bloch oscillation cycle by getting 
its internal state rotated somehow. In the perspective 
of the cat state, a cycle of Bloch oscillation is a cycle 
of birth-growth-death. The cycle can be interrupted by 
putting _F = at an appropriate time, e.g., at t = T/2, 
when the cat is in its largest size. Or if the direction of 
the force F is reversed at i = T/2, the two wave pack- 
ets will continue moving rightward instead of going back. 
This would help to increase the size of the cat further. 

We have followed the center-of-mass motion of the 
wave packets in time. The results are shown in Fig. [3] 
(solid lines). We see that the semi-classical theory (dot- 
ted lines) is correct quantitatively. The numerically exact 
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FIG. 3: (Color online) Evolution of the center-of-mass of the 
wave packets \['i,o(a;,i) in a Bloch oscillation cycle. The pa- 
rameters are the same as in Fig. [T] and Fig. [21 The results 
obtained by solving Eq. ([5]) numerically/exactly (solid lines) 
agree very well with those (dotted lines) by the semi-classical 
theory [see Eq. ([J}]. 



results deviate significantly from the semi-classical pre- 
dictions only in the vicinities oi t = and t = T. The 
reason is that due to the sudden change of 0, both wave 
packets are partially excited and some packets (not visi- 
ble in the snapshots) belonging to higher bands are emit- 
ted around t = 0, which may perform Bloch oscillation 
also and return around t — T . 

We now turn to the problem of the feasibility of the 



scheme in experiment. In our simulation, Fa/Er = 
0.005. It is essential to make sure that this ratio is much 
smaller than unity. First, the potential drop Fa between 
two neighboring sites should be much smaller than the 
gap between the zeroth and first bands (see Fig. [T]d) so as 
to suppress Zener tunneling [2[ . Or equivalently, the Bril- 
louin zone, especially its boundary, should be traversed 
slowly so that transition into excited bands can be ne- 
glected. Second, according to (O, the distance between 
the two wave packets is inversely proportional to F, thus 
smaller F means larger separation or larger "cat". Of 
course, there should be an optimal value of F since the 
period T is also inversely proportional to F. For ^^Rb 
and a lattice constant a = A/2 = 392.5 nm, the ratio 
above corresponds to a tilt angle ip — 4° and a period 
T = 53 ms. On the contrary, under the chosen detun- 
ing and strength of the optical lattice, the spontaneous 
radiation rate T^ff of the atom is about 0.2 s~^. Thus 
the atom is long lived enough to oscillate several cycles 
before incoherent processes set in. 

In conclusion, we have proposed that the spin- 
dependent optical lattice may offer an opportunity to cre- 



ate Schrodinger cat states by using Bloch oscillation. Our 
scheme has several interesting advantages. First, the cat 
state experiences birth-growth-death cycles repeatedly. 
It would be worthy to study experimentally how this pro- 
cess is damped in a real optical lattice, which is believed 
to be well isolated from the environment. Second, if we 
start from a Bose-Einstein condensate and minimize the 
atom-atom interaction which is deleterious to the Bloch 
oscillation, it might be possible to create a collection of 
atoms condensed in a cat state. We note that some gen- 
eralizations are also possible. For example, though here 
we focused on the one dimensional case, the scheme can 
be directly extended to higher dimensions [l^ since two- 
dimensional spin-dependent optical lattices have already 
been demonstrated experimentally |29| . Furthermore, in 
contrast to the static force considered here, periodically 
modulations [la . Il7| are worth consideration also since 
they may help to increase the size of the "cat" . 
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